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Abstract 

We use a simple model to describe the nonlinear dynamics of a dense two dimensional dipolar 
exciton gas. The model predicts an initial fast expansion due to dipole-dipole pressure, followed 
by a much slower diffusion. The model is in very good agreement with recent experimental results. 
We show that the dipole pressure induced expansion strongly constrains the time available for 
achieving and observing Bose-Einstein quantum statistical effects, indicating a need for spatial 
exciton traps. We also suggest that nonlinear ballistic exciton transport due to the strong internal 
dipole pressure is readily achievable. 
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A new species of exciton, the spatially indirect dipolar exciton, has recently emerged in 



the continuing search for degenerate quantum gases and liquids in excitonic systems 
Spatially indirect dipolar excitons are coulomb bound electron-hole pairs in which the elec- 
tron and the hole are confined in two spatially separated two-dimensional (2D) layers. Such 
excitonic states are usually achieved in double quantum well (DQW) structures with an 
external electrical bias apphed perpendicular to the quantum well (QW) plane to lift the 
degeneracy of the confined states in the two respective QW's . As a result, the lowest 
energy optical transitions occur between the two QW's, leading to the formation of indirect 
excitons. While the binding energy of such an exciton is only slightly modified due to the 
spatial separation, the lifetime is enhanced by many orders of magnitude compared to di- 
rect excitons (where the electron and the hole are in the same layer). This extremely long 
lifetime makes them especially interesting in the context of degenerate quantum gases of 
excitons, since it may provide sufficient time for thermalization as well as cooling of initially 
hot excitons. Another prominent feature of this type of indirect exciton is its permanent 
dipole moment, aligned perpendicular to the QW plane, due to the spatial separation of 
the electron and the hole. This dipole moment gives rise to a repulsive dipole-dipole force 
between excitons. On one hand, this repulsion prevents the formation of e-h complexes such 
as biexcitons and e-h droplets, and therefore makes such systems even more attractive can- 
didates for observing degenerate Bose- Einstein gases. On the other hand, it leads to a strong 
driven expansion of a dense dipolar exciton gas, as was very 

this dipole-dipole repulsion seems to constrain the time that such an indirect exciton gas 
can remain dense after photoexcitation, an obvious obstacle to the observation of quantum 
statistical effects of excitons. 

Here, we use a simple model to describe the dynamics of such dipolar excitons, taking into 
account the dipole-dipole interaction. Besides quantitatively describing current experimen- 
tal observations of the nonlinear expansion dynamics of dipolar exciton gases, this model 
suggests that the expansion caused by the dipole repulsion indeed strongly limits the time 
that a high mobility exciton gas remains dense, explaining the difficulties in the observa- 
tions of exciton condensation. Some sort of spatial confining (exciton traps) might then be 
necessary for achieving a dense and cold dipolar exciton gas, a prerequisite for excitonic 
Bose-Einstein condensation in such systems. 

Two recent experimental observations of dipolar exciton dynamics by Voros et al. 
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(Ref. [3|) are the basis of our model assumptions: (a) The dipolar exciton expansion dy- 
namics can be divided into two regimes: immediately after the laser excitation pulse, a very 
fast expansion is observed, in which the excited exciton cloud expands from its initial small 
diameter (~ SOfim) by a factor of ~ 6, in a few tens of nanoseconds. At later times, the ex- 
pansion slows down significantly, and can be very well described by exciton diffusion process, 
(b) The diffusion coefficient of the dipolar excitons, Dx, has a power law dependence on the 
well width, L, where it is found to be of the form Dx oc L^. This specific dependence was 
found previously for free electrons in similar structures by Sakaki et al. in Ref. 141, which 
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is well explained by interface roughness scattering from well width fluctuations 
This finding is a strong evidence that for dipolar excitons at low temperatures, potential 
fluctuations due to interface roughness is also the dominant scattering mechanism. 

Based on these two observations, we suggest a simple model to describe the dynamics 
of the dipolar excitons (see also Ref. [8|). The model includes the dipole-dipole repulsion 
between the excitons as the main driving force for the initial fast expansion, and the static 
interface roughness scattering as the drag mechanism. 

The dipole-dipole interaction energy per exciton of a gas of perfectly aligned dipolar 



excitons is given by 



£dd{r) = — - — nx = anx, (1) 



where nx = nx{r) is the 2-dimensional density of the dipolar excitons, Zq is their dipole 
length, and e is the background dielectric constant. The resulting dipole-dipole repulsion 
force, 

Fdd = -aVnx, (2) 

proportional to the density gradient, will drive the exciton expansion. The other effective 
force that is responsible for the diffusion of an exciton gas is 

Fd = -VC- (3) 

Here ( = kTxln (l — e~^°/-^-^) is the chemical potential in the non-interacting limit j^, where 
To = is the degeneracy temperature and Tx and mx are the exciton gas temperature 

and exciton effective mass, respectively. The overall internal force driving the excitons is 
therefore F = Fad + Fo- The resulting exciton current J = rix/iF = nxfiiFdd + FD) consists 
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of contributions from both the dipole-dipole repulsion and the diffusion: 



Jd = -nxfJ^aVnx (4) 

kT 

Here, /z is the exciton mobihty which is, in general, a function of |F| = F. 

We now construct a simplified model to obtain the dependence of fi on F following 
the observation that interface roughness scattering is the dominant relaxation mechanism. 
We model the interface roughness scattering by a random distribution of quasi-elastic, hard 
sphere scattering centers with a characteristic mean free path I. In the low-field limit, the mo- 
bility, fiif, is independent of F. In that limit, fiif is related to the diffusion coefficient through 
the modified Einstein relation^] D = iikTx / {e^°^^^ — 1), which in the non-degenerate limit 
reduces to the well-known jj^if = Dx/kTx- In this picture, Dx = l"^ /t = Ivth, where r is 



the exciton scattering time and Vth is the thermal velocity Vth = ^/'^kTx/rax- In the high 
field limit however, the drift velocity Vd depends on the force F and is determined by the 



average velocity the exciton gains between successive collisions, given by Vd = a/ Fl/{2mx)- 
The resulting high field limit for the mobility is then fihf = Vd/F = \Jl/ {2mxF). We thus 
construct a smooth function for the mobility as a function of F that conforms to the above 
two limits: 

The crossover between the two transport regimes occurs at an F which corresponds to 
Vth ~ 4frf. The dependence of the mobility on Vd/vth is plotted in Fig. [1^. The free exciton 
gas dynamics is now described by a nonlinear diffusion-drift equation, 

^ + V • (J,, + J,) + ^ - /x(r,t) = 0, (7) 
at Tx 

where Ix{f-,t) is the exciton source and tx is the exciton recombination lifetime. This 
nonlinear equation can be solved numerically. 

Since /i depends on the internal force F, the exciton transport is in general nonlinear, 
analogous to nonlinear I-V characteristics of charged particles. In the case that Vd <^ Vth, 
H becomes independent of F and the exciton transport becomes linear. However, even in 
the linear transport limit, the exciton expansion can be very fast, strongly driven by the 
dipole repulsion which can dominate over diffusion (it turns out that this is the case in the 
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experiments presented in Ref. j3[). Thus, we define the following quantity to describe the 
importance of the fast dipole driven expansion: 

7(r, t) = J J Jo = ^, (8) 

which, in the dilute limit, reduces to just the ratio between the dipole-dipole energy term 
and the thermal energy of the excitons, 7 = When 7 3> 1, the expansion is driven by 
the repulsive dipole-dipole force. When 7 <^ 1, the expansion becomes diffusive. An exciton 
gas will exhibit a fast driven expansion if its initial density profile, = nx{r = 0,t = 0), 
satisfies 7(r = 0, t = 0) > 1. As a dense exciton gas starts to rapidly expand due the dipole 
forces, its density will decrease and the expansion will slow down. Wherever the density 
decreases such that 7(r = 0, t = 0) < 1, the expansion will become diffusive, until eventually 
the diffusion will dominate the dynamics of the whole exciton gas. 

We start our discussions by comparing the model calculations to the experimental results 
of Ref. ^. There, DQW structures were used with different well widths: 80A, lOOA, 120A, 
and 140A and the corresponding Dx were found to be (in units of cm'^/s): 0.24, 0.74, 2.08, 
and 9.4. These numbers yield a corresponding mean free path of l{nm) = 1.4, 4.5, 12.5, 56.6, 
for the experimental temperature, T = 2K, and for an exciton effective mass rrix = 0.2me. 
Note that for the wide DQW structures I > ax (where ax ~ lAnm is the exciton Bohr 
radius) while for the narrow DQW's, I < ax- Using the parameters of Ref. Q], we find 
that exciton drift velocity Vd is always smaller than Vth even for samples with the highest 
Dx- The exciton transport is therefore linear with /i ~ However, we estimate that 
7(r = 0, t = 0) ~ 800, indicating that the exciton cloud should exhibits an initial fast driven 
expansion due to the dipole pressure, and later a diffusive behavior. 

For our calculations, we assume radial symmetry and an initial gaussian distribution of 
optically excited dipolar excitons, with a half width ctq = 15fim and an = lO^^cm"^ 
(unless stated otherwise), similar to the experimental conditions of Ref. We also use an 
exciton dipole length of zq = 12nm. A temperature of 3K and an exciton lifetime of are 
assumed. Fig. ^ presents calculated exciton density profiles for several different times and 
for Dx = 9Acm'^/s, showing the exciton expansion. 

Fig. |2k plots the calculated squares of the half-width - half-maximum (a^) of the exciton 
density profiles at various times, for various Dx extracted from the experiments of Ref. P| . 
Both the fast initial expansion and the later slow diffusive regime are clearly seen. The 



larger the Dx, the faster the initial expansion, as /i and Dx are linked through the scattering 
length. While we did not try to do any "fine tuning" of parameter fitting, there is a very 
good agreement with the experimental results (Fig. 4 in Ref. The only parameters for 
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FIG. 1: (a) Normalized dipolar exciton mobility as a function of the ratio of the drift velocity to 
the thermal velocity Vd/vth- (b) Dipolar exciton radial density distribution profiles for different 
times. 

The experimental data, however, seem to have a sharper transition from the driven regime 
to the diffusive one than the calculation. This is seen in the sharper bend of cr^ with time 
shown in Ref. [^. This discrepancy may be due to the heating of the exciton gas as it 
expands, transferring the potential energy into thermal energy by the elastic scattering 
processes. Such heating process is not included in our model. 

Fig. 1213 presents the calculated time traces of the emission from different fixed radial 
positions of the expanding exciton gas (for Dx = 0.74cm^/s). These calculated time traces 
are also in a good agreement with the experiments (Fig. 2 in Ref. jsj). It seems then that 
this model can quite accurately approximate real experimental observations, based on the 
few simple assumptions given above. 

After establishing the validity of this picture, we are now ready to discuss some of its 
important consequences. It is already clear that a dense and mobile dipolar gas will tend to 
expand quickly due to the internal repulsive forces. If the initial cloud is small, the expansion 
will quickly reduce the exciton gas density and will limit the time the high density condition 
can be sustained. We start by discussing small initial exciton clouds, excited by tightly 
focused laser beams, which is the most common experimental approach. In Fig. we plot 
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FIG. 2: (a) cr^ of the density distribution profiles as a function of time, for the diffusion coefficients 

of Ref. (b) Calculated dipolar exciton emission density as a function of time from different 

radial positions of the exciton cloud, corresponding to the experimental positions measured in 
n 

Ref. jj]. The parameters for the calculations in both (a) and (b) are given in the text. 

the central density of an expanding, initially small (cr = 15/im), high mobility exciton gas 
{Dx = lOcm^/s). The calculations are performed with a gas temperature of Tx = 2K, as 
a function of time, for different initial exciton densities n^. the horizontal dash-dotted line 
represents the density, n^^, in which the exciton thermal de-Broglie wavelength is equal to 
the exciton inter-particle distance, i.e., ^^dB^'x/dx = 1- One expects to observe quantum 
statistical effects due to Bose-Einstein distribution only for densities larger than Ux- As can 
be seen, within less than 30ns, the density drops below Ux even for high initial gas densities. 
This time is much shorter than the exciton lifetime of 3/is. Again, this is a consequence of 
the fast initial expansion. Fig. [Hb shows the similar density-time traces for an initial density 
= lO^^cm"^ for different temperatures. The circles on each density-time trace mark Ux 
for the particular temperature, given by: rix = "'T^'XgkTx/i^Trh'^), and the curved dashed-dot 
line is a guide to the eye. While decreasing temperature results in an increased mobility 
for a given Dx, there is also a linear decrease in rix, resulting in an improvement of the 
time interval in which nx > n\. However, even at Tx = 0.57^ this time is still under 70ns, 
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significantly shorter than the exciton radiative hfetime, tx- Fig. I^fc shows similar traces for 
different diffusion coefficients. A smaller diffusion coefficient is related to a smaller mean 
free path /, which in turn yields less expansion and less density decrease. 
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FIG. 3: Center exciton densities as a function of time, nx{t,r = 0) for (a) different initial center 
densities n^, (b) different exciton gas temperatures Tx, and (c) Different exciton diffusion coeffi- 
cients Dx- The relevant parameters for the calculations in (a), (b) and (c) are given in the text. 
The dashed-dot lines in (a) and (c) and the circles in (b) mark the de-Broglie density n'^. 

Fig. |3^,b summarize the above results. The time interval At in which nx > n-x as a 
function of the exciton gas temperature, for various Dx values is shown in Fig. and for 
various initial densities Ux in Fig. |3Jd. For high mobility (diffusivity) samples, decreasing 
the temperature and increasing the initial density only slightly increases At, to values still 
much smaller than tx- This is due to the fast driven expansion. For low mobility samples, 
at very low temperatures. At can be in the microsecond range, comparable to tx- For these 
samples, however, the mean free path, /, is much smaller than the exciton Bohr radius, ax- 
This may complicate the observation of Bose-Einstein quantum statistical effects in general 
and excitonic BEC in particular due to the disorder induced localization of the excitons. 
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FIG. 4: Time interval At in which the center dipolar exciton gas density, nxir = 0) is larger than 
the de-Broglie density, function of the exciton gas temperature, Tx- (a) shows the results 

for different diffusion coefficients, Dx while (b) shows different initial center densities, n^^. (c) 
shows A.t as a function of the initial size of the exciton cloud for Tx = 2i^. 



Fig. Ht shows the dependence of At on the initial size of the exciton cloud or Dx = 
lOcm^/s and Tx = 2K. As the initial size increases (keeping the same center density), the 
dipole driving force decreases Fdd oc |Vrax| oc I/ctq, and so does the expansion rate. Starting 
with large enough initial exciton clouds leads to a much smaller expansion rates and thus 
to At values comparable to the exciton lifetime tx- 

To conclude this analysis, there are two strategies to keep a high mobility dipolar ex- 
citon gas dense enough for long enough time intervals. The first is to optically excite an 
initially very large cloud, thus suppressing the fast expansion. This approach, however, is 
not favorable since it involves high excitation power which could lead to complications such 
as increased heating of the sample. The other strategy is to design an artificial spatial con- 
finement in the QW plane. For such dipolar excitons, the trapping potential can be induced 
«the. through a local strain m the QW Q or by locally applymg the external electrostatrc 
potential[ll]. These exciton traps not only keep the excitons dense, but also eliminate the 
extra heating (due to the dipole-dipole potential) of a free-expanding exciton gas mentioned 
previously and breaking the 2D translational which allows for the the formation of BEC. 

It is interesting to compare the two characteristic densities of this problem. The first is 
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the de-Broglie density, that marks the onset of a distinguishable Bose- Einstein statistics 
effects. The second density is the critical density for the driven expansion regime, = 
kTx/o., derived by setting 7 = 1 in the linear transport limit and assuming non-degenerate 
gas. ra^ marks the transition from a fast dipole-pressure driven exciton expansion to a 
diffusive expansion. Their ratio yields, 

n'x _ 2mxge^ZQ 



rf^ he 



20, (9) 



for typical GaAs DQW structures and is independent of temperature, justifying the above 
assumption. This means that a degenerate Bose-Einstein gas of dipolar excitons is always 
strongly driven by the dipolar pressure. This can have some interesting consequences. For 
example, it indicates that a highly degenerate dipolar gas will expand quickly, in contrast to 
most atomic BEG gases that tend to expand very slowly. Also, it is not clear how a driven 
expansion of a highly degenerate exciton gas would differ from that of a classical gas and 
how would the scattering of such an exciton gas from well width fluctuations would differ 
from the simplistic, ballistic "billiard ball" scattering model we have used. Our model also 
suggests that for high mobility samples, such as the ones used in Ref. Q], the mobility is on 
the edge of the nonlinear transport regime. By cooling the exciton gas further or by growing 
wider samples the predicted range for the study of highly nonlinear transport due to the 
internal dipole pressure can be reached for the first time. 

We thank Loren Pfeiffer, P.M. Platzman, David Snoke and Xing Wei for enlightening 
discussions. 



[1] D. Snoke, Science 298, 1368 (2002). 

[2] L. V. Butov, J. Phys.: Condens. Matter. 16 R1577 (2004). 

[3] Z. Voros, R. Balili, D.W. Snoke, L. N. Pfeiffer, and K. West, Phys. Rev. Lett. 94, 226401 
(2005), |cmId-mat/050415T] (2005). 



[4] H. Sakaki, T. Noda, K. Hirakawa, M. Tanaka, and T. Matsusue, Appl. Phys. Lett. 51, 1934 
(1987). 

[5] R.E. Prange and T. W. Nee, Phys. Rev. 168, 779 (1968). 
[6] A. Gold, Solid State Commun. 60, 531 (1986). 
[7] A. Gold, Phy. Rev. B 35, 723 (1987). 

10 



[8] A. Ivanov, P.B. Littlewood and H. Haug, Phys. Rev. B 59, 5032 (1999); A. Ivanov, Europhys. 

Lett. 59(4), 586 (2002). 
[9] S. Ben-Tabou de-Leon and B. Laikhtman, Phys. Rev. B 63, 125306 (2001). 
[10] D. Snoke, Y. Liu, Z. Voros, L. N. Pfeiffer, and K. West, Solid State Commun. 134, 37 (2005). 
[11] R. Rapaport, G. Chen, S. Simon, O. Mitrofanov, L. N. Pfeiffer, and P. M. Platzman, Phys. 
Rev. B, 72, 000000 (2005), |cmldiai:t70504178| (2005). 



11 



